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Abstract 

We explore the tunneling dynamics of strongly correlated bosonic mixtures in a one-dimensional 
double-well. The role of the inter- and intra-species interactions and their interplay is investigated 
using the numerically exact Multi-Configuration Time Dependent Hartree (MCTDH) method. The 
dynamics is studied for three initial configurations: complete and partial population imbalance and 
a species separated state. Increasing the inter-species interaction leads to a strong increase of the 
tunneling time period analogous to the quantum self-trapping for condensates. The intra-species 
repulsion can suppress or enhance the tunneling period depending on the strength of the inter- 
species correlations as well as the initial configuration. Completely correlated tunneling between the 
two species and within the same species as well as mechanisms of species separation and counterflow 
are revealed. These effects are explained by studying the few-body energy spectra as well as the 
properties of the contributing stationary states. 
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I. INTRODUCTION 



Ultracold atoms represents a distinct phase of matter for the exploration of fundamental 
quantum processes {l-3|. The ability to precisely control ultracold systems has triggered 
investigations in different fields such as quantum simulation and information processing 
4|, quantum phase transitions 5, 6| and driven quantum systems Experimentally it is 
possible to control not only the external potentials but also the effective interactions between 
the atoms using Feshbach resonances (Sj. Moreover, dimensionality can be tuned and in 
particular quasi one-dimensional systems can be achieved by confining two transverse degrees 
of freedom. In such waveguide-like systems, confinement induced resonances ^| provide an 
additional tool to tune the interactions thus making the study of strongly correlated systems 
experimentally feasible. An interesting observation for a single species bosonic system in one 
dimension, is that infinitely strongly repulsively interacting bosons possess the same local 
properties as a system of non-interacting fermions. This effect, known as fermionization has 
3een experimentally observed 26|, |27| and can be explained via the Bose-Fermi mapping 



28| 



Inspired 
perimental 



by the results of single bosonic species, recently, there has been a lot of ex- 



1Ch14| and theoretical 15N25| interest in the static properties of multi-species 



bosonic mixtures. In these systems, the interplay between the inter- and intra-species forces 
as well as different masses or potential asymmetry gives rise to various phenomena and ef- 
fects not accessible in the single component case. For instance, the process of composite 
fermionization occurs when the inter-species coupling is set to infinity and the strong repul 



sion provides different pathways for phase separation [161, [171, |21|. Moreover instabilities [l5| 
as well as new phases such as paired and counterfiow superfluidity 25| have been observed. 

Focusing on the quantum dynamics, the double well provides the simplest prototype for a 
finite lattice and Josephson junction and is especially a very elucidating case for studying the 
fundamental characteristics of quantum tunneling. Theoretically, the tunneling dynamics of 
single species through the crossover from weak to strong interaction regimes reveal interesting 
effects such as Josephson oscillations, pair tunneling, self trap ping as well as fermionization 
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46| which have also been observed experimentally 32l. |33|. 



More recently these studies have been extended also to systems of binary bosonic mixtures 



34|-|4l|. These works demonstrate various effects such as macroscopic quantum self-trapping 



namics 
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and coherent quantum tunneling 3J|, observations of collapse and revival of population dy- 



40] . symmetry breaking and restoring scenarios 37] as well as dipole oscillations 



induced pairing and counterflow superfluidity 



4l|. 



However, most of the work has been done on the mean-field level either by solving Gross- 
Pitaevskii equations or by using the lowest band Bose-Hubbard model. Although these 
studies do provide interesting insights into the mechanism of tunneling, an investigation of 
the complete crossover from the weak to strong interaction regime allows the examination 
of new effects and mechanisms not present e.g. in the mean field description. For instance 
referring to the case of two species in a harmonic trap, it has been found that if one species is 
localized due to its heavy mass then it can act as an effective material barrier through which 



the lighter component tunnels 
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43j. The feedback of this material barrier leads to different 



pairing mechanisms for the light species. Moreover, few body systems provide a bottom-up 
approach towards the understanding of many-body phenomena. Experiments exploring few 



atom systems in finite optical lattices 



44] serve as promising setups for designing transistor- 



like structures from the perspective of atomtronics. 

In this paper we study the tunneling dynamics of a binary mixture of bosonic species in 
a one-dimensional double-well from a few body perspective. Using the numerically exact 
Multi- Configuration Time-Dependent Hartree method (MCTDH, see Appendix E]) 471. l48] . 
we investigate the crossover from weak to strong interactions focusing in particular, on 
microscopic quantum effects and mechanisms which are prominent in few body systems. We 
demonstrate how the interplay between the inter- and intra- species interactions affect the 
rate and behavior of the tunneling in a non-trivial way. A strong increase of the tunneling 
period is observed as the inter-species repulsion is increased. However in certain cases, 
and especially for the strongly-interacting regime, increasing the intra-species interactions 
leads to an increase of the tunneling rates in contrast to what is observed for single species 
systems. Preparing different initial states leads consequently to a diverse tunneling behavior. 
For complete imbalance of the populations, i.e., when the particles are all prepared initially 
in the same well, or when the species are localized at different wells (species separation) the 
tunneling is strongly correlated meaning that the species tunnel either in phase or out of 
phase . Only for very strong intra-species interactions these correlations are reduced. On the 
other hand for partial population imbalance e.g. one species is delocalized and the second 
one is localized, a mechanism of species separation and counterflow appears. The various 



effects are attributed to the features of the energy spectrum and explained by examining 
the density profile of the contributing stationary states. 

The paper is organized as follows. In Section [Til we introduce our model and setup. 
Subsequently we present and discuss the results for the quantum dynamics of the mixture 
with three bosons (two bosons of species A and one of species B). Three initial state scenarios 
are examined: complete population imbalance in Sec. IIIH complete species-separation in 
Sec. IIV| and partial imbalance in Sec. |Vl The computational method MCTDH is described 
in the Appendix Rl 



II. MODEL AND SETUP 



We consider a mixture of two species of bosons labeled by A and i? in a one-dimensional 
double well potential. These may correspond to two different kinds of atoms or could be two 
hyperfine states of the same atomic species. The fact that there are two different species 
induces distinguishability and thus fundamentally alters the physics and in particular the 
quantum dynamics compared to the case of a single species. 

Our Hamiltonian reads (see j21| for details) 



(T=A,B i=l 



Pl,i 



i<j 



Na Nb 

^^^AB(xA,i -a:B,i)- (1) 
i=i j=i 



where Ma,b is the mass for species A and B, respectively. 

We assume here that the different species obey the same single particle Hamiltonian, 
i.e., they possess the same mass and experience the same single-particle potential. The 
double-well potential U{x) = ^Mux"^ + h6co{x) is modeled as a harmonic potential with a 
central barrier shaped as a Gaussian h6^{x) = h '^ of width s = 0.5 and height h = 8.0. 

Dimensionless harmonic-oscillator units i.e.. Ma = Mb = 1, u = 1 are employed throughout. 
In the ultracold scattering limit, one can approximate the interaction (both intra- and 
inter-species Vab) with an effective contact potential 



Numerically we sample the delta-function as a very narrow Gaussian (choosing of course a 
spatial grid dense enough to sample this narrow peak). 

The different initial configurations are achieved by adding a tilt to the double-well which 
can be different for the two species depending on the required state. Thus an individual 
well could be made energetically more favorable (tilted) for a certain species. For instance, 
to prepare a complete imbalance, the double wells of both species are tilted the same way, 
while to prepare a species-separated scenario Ua is tilted opposite to Ub- To prepare the 
partial population imbalanced state one has to tune the tilt for both species judiciously 
depending on the given interaction strength such that the required population configuration 
is achieved. The ground-state is then computed by the relaxation method and results in the 
desired initial state. For the study of the dynamics the tilt is instantaneously ramped down 
to obtain a symmetric double-well at t = 0. 

In order to investigate systematically and in detail the tunneling processes for binary 
mixtures we consider the simplest non-trivial few-body system consisting of two bosons of 
species A and one of species B. This system captures the most important microscopic 
quantum dynamical processes occurring for few-body bosonic mixtures. In this case we 
have two independent parameters qab and qa (since there is only a single boson B species). 
When the inter-species interaction qab is zero, the two components are completely decoupled 
meaning that the single B boson will undergo Rabi oscillations between the wells. The A 
bosons will then follow a correlated two-particle dynamics regulated by the intra-species 



interaction q 



literature 



45 



?his case is not addressed here but has been discussed in detail in the 
Another case which reduces to that of a single species is qab — ^ gA, 
where the essentials of the tunneling dynamics is that of three particles of a single species. 
Our focus is exclusively onto the cases where we expect significant deviations from the single 
species scenario. 



III. COMPLETE POPULATION IMBALANCE. 

We begin our study by exploring the quantum dynamics for an initial state where all the 
atoms are loaded into the left well. As observables, we compute the time evolution of the 
one-particle density of each species and the resulting population in each well. For the right 
well we have 
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Table I. Tunneling periods for different values for the case of complete imbalance. 



Period 


QAB = 0.0 


QAB = 0.2 


QAB = 5.0 


QAB = 25.0 


9A = 0.0 


2 X 10^ 


1 X 10^ 


9 X 10^ 


1 X 10^ 


9A = 0.2 


6 X 10^ 


6 X 10^ 


1 X 10'' 


1 X 10*^ 



poo 

na{t) = Na pa{x;t)dx (2) 
Jo 

where is the one-body density of the species a = A,B and the total population of the 
right well is ur = ha + ns- Due to symmetry and resonant mechanisms, we always have a 
complete transfer of the population of both species between the two wells which happens in 
most cases according to a periodic pattern with period T. 



1. Repulsive interspecies interaction and binding mechanisms 

The most important effect of increasing the interspecies interaction qab is a very strong 
increase of the tunneling period up to very large values. This can be seen in Table HI where we 
show the tunneling period with increasing Qab- This behavior is counter intuitive since with 
increasing repulsion between the species initially localized in the same well, one would expect 
the tunneling to be enhanced. The delayed tunneling is reminiscent of the one found for the 
case of a single species j2l[ I^gJ and is the few-body equivalent of self-trapping. The primary 
reason for this decrease of the tunneling frequency, especially for low interactions (within 
the so-called Bose-Hubbard regime) can be attributed to the energy spectrum presented in 
Fig. [1] considering the states that contribute to the dynamics. The eigenstates are typically 
characterized by the superpositions of different number states such as \AA,B), where the 
vector indicates two A boson occupying the left well and one B boson in the right well. As 
qab increases, different doublets are formed in this lowest band. The energetically highest 
doublet shown in Fig. [H consisting primarily of the states \AAB^{]) ± \{)^AAB), is of 
relevance to our case since these eigenstates possess maximum overlap with our initial state 
\AAB,Q). 

With increasing Qab-, the number states |y4y4i?, 0) and |0, AAB) depart energetically from 
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Figure 1. (color online) Energy spectrum for qa = 0.0 and small interaction strengths gAB- 



other number states due to their big on-site interaction energy, (having all bosons in the 
same well) and the eigenstates \AAB, 0) ± |0, AAB) become increasingly degenerate thereby 
forming a doublet. The tunneling then consists oscillations between | AAi?, 0) and |0, AAB), 
while the decreasing energy splitting of the doublet leads to an increase of the tunneling 
period. This is the few-body analogue of the self trapping mechanism in single-species 
condensates. The impressive fact though, is that this behavior is even more pronounced for 
higher interactions (see eg. qab = 5.0, 25.0 in Table |T]), in contrast to the single species case 



fsee ref 



21 



46|) where there is a reduction of the period due to higher band contributions 
and fermionization (with the particles tunneling as uncorrelated fermionized bosons with 
the corresponding Rabi frequency) in the strong interaction limit. In our case, as long as 
the interaction takes place predominantly between the different species (here for simplicity 
Qa = 0), firstly there is no fermionization in the regular sense, and secondly, the particles 
tunnel in a highly correlated manner - meaning that the initial localized state \AAB,0) 
does not tunnel to states like \B, AA) or \AB, A) since they possess much lower interaction 
energy and are thus energetically off resonant. The fact that we encounter correlated particle 
tunneling (all bosons together) is documented by the probability of finding all the particles 
in the same well which remains very close to unity throughout the dynamics. As a first 
conclusion, we see that the repulsive inter-species interaction causes 'binding' between the 
particles and reduces the tunneling rates. 
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Figure 2. (color online) Population of the right well (a) ua of species A and (b) of species B at 
9AB = 0.2 for different values. 

2. Intra-species interactions to control tunneling and correlations 

In contrast to the one-way effect of the qab in the tunnehng dynamics, the intra-species 
repulsion qa plays a more complicated role - controlling both the tunneling period and 
the degree of correlations. Let us first explore a weak inter-species interaction strength 
Qab = 0.2. In Fig. |21 we illustrate the tunneling dynamics for different values of qa at 
9 AS = 0.2 for species A and B by showing the population of the right well nA,nB- Except for 
the periods T of the tunneling envelope, one also observes rapid small amplitude oscillations. 
Concerning T, we obtain a monotonic increase as qa = 0.0 0.3. However this behavior 
changes as we go beyond the weak interaction regime for qa and we observe a decrease of 
the tunneling period for qa = 5. Another important feature is that the two components A 
and B undergo roughly the same evolution of the oscillation pattern (compare Fig. [2] (a) 
and (b)) which is suggestive of strong inter and intra-species correlations in the sense that 
all bosons tunnel together. This changes slightly only for very strong interaction qa = 25, 
where the tunneling period reduces substantially while the pattern becomes more erratic 
consisting of two primary oscillations and unlike the previous cases, the dynamics of the two 
components is not completely identical. This indicates, in the line of argumentation provided 
above, a reduction of the correlations between the two species and attempted single-particle 
tunneling. 

Within the weak interaction regime where the effective lowest band number states de- 
scription is valid, the tunneling process of shuffling between the two completely local- 
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ized states \AAB,0) and \0,AAB) can be more specifically described by the sequence 
\AAB,0) \AB,A) \B,AA) \0,AAB). The effective tunneling rate within this 
lowest band description, is in general approximately given by 

f j'/{Ei-E2){Ei-Es) (3) 

where J is the effective coupling term between the two sites and Ei is the energy of the 
initial and final number-state of the tunneling and E2, E3 are the energies of the intermediate 
number-states respectively, which in this case are \AB,A) and \B,AA). Considering the 
interaction part of the number states, the completely localized states 0) and |0, AAB) 

have energies ~ 2gAB + gA while the states \AB,A) and \B,AA) have energies ~ qab and 
qa respectively. Therefore, for this tunneling process, the tunneling rate according to Eq. 
([3]) scales as / ~ I'^gABidAB + 9a)- From this relation follows that the tunneling rate 
decreases for increasing qa in the weak interaction regime as we have seen above. 




Time Time 

Figure 3. (color online) Population in the right well (a) ua of species A and (b) of species B at 
QAB = 5.0 for different qa values. 

The opposite effect i.e. decrease of the tunneling period, as qa increases further beyond 
the weak coupling regime, is attributed to the increasing splitting of the main contribut- 
ing doublet \AAB,Q) ± \{],AAB) as qa increases. Additionally, higher bands contributions 
appear, especially for stronger intra-species interactions qa = 25.0, which breaks the com- 
pletely correlated tunneling behavior, allowing for attempted single particle tunneling into 
energetically higher number-states like \AB^A). 

Turning now to higher inter-species interaction qab = 5.0 we observe in Fig. [3] that the 
tunneling period decreases strongly as qa increases. Since qab is in this case beyond the 
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Figure 4. (color online) One-particle density as obtained by taking into account the most important 
eigenstates contributing to the initial state for (a) species A and (b) species B for qab = 5.0. 

weak coupling regime, we focus on an analysis of the density profiles of the contributing 
eigenstates shown in Fig. IHto understand the effect of increasing qa- As qa increases, the 
repulsion of the A bosons leads to a broadening of their density profile. This broadening 
leads to a greater overlap of the wave functions of A atoms localized in the left and the 
right well and this in turn increases the effective tunneling coupling and the corresponding 
tunneling rates. At qa = 0, the localized densities pA and ps are spatially separated in 
each well as a consequence of the repulsion between the species. Note that the density of 
the B boson possesses its maximum for larger values of |x| thereby 'sandwiching' the A 
boson population. This arises from the fact that due to the unequal number (A^^ > Nb), 
it is energetically favorable to shift the density of the B species to larger values of 
As Qa is increased, the two localized densities Pa,Pb in the two wells gain an increasing 
overlap which can be observed as a vertical upward shift of the density profile at x = that 
becomes progressively stronger with increasing qa- This mechanism, also present for other 
contributing states, leads to an overall increase of the tunneling coupling and consequently 
to an increase of the tunneling frequency for strong interactions. 

The overall features with respect to the different time-scales and oscillatory tunneling 
behavior is similar for very strong interspecies interactions qab = 25.0 with the exception 
of Qa = 20.0 [Figj5]^a)]. Only this case can be considered as a tunneling mechanism close to 
fermionization. In this regime, the bosons become isomorphic with non-interacting fermions 
and thus the tunneling dynamics approaches that of independent non-interacting fermions. 
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Figure 5. (color online) (a) Population in the right well of the total hr and the individual species 
UA, ns for gA_B = 25.0 = 20.0. (b) Pair/triple probability. Paa, Pab and Paab correspond to 
the probability of finding AA, AB and AAB boson in the same well, respectively. 

The latter tunneling frequency is close to the Rabi-frequency which is significantly faster than 
the previously discussed cases. To understand the reduction of period from a number-state 
perspective we note that the nearly fermionized bosons occupy both the lowest band and 
the first excited band. As a consequence, the previously off-resonant intermediate number- 
states (namely the states |AAi?,0), \AB,A) and \B^AA)) become near resonant, since the 
particles can tunnel between the excited band of the two wells without significant change of 
energy. This results in a reduction of the effective tunneling period while the separation of 
the time-scales involved in the dynamics is strongly reduced. The correlations among the 
bosons with respect to the tunneling process is also strongly reduced as can be seen in Fig. 
|5]^b), where the pair and triple correlations show a strong deviation from the value 1, i.e. 
from the strongly correlated case. 



IV. SPECIES-SEPARATED INITIAL STATE. 



Let us now consider the initial state for which the two species are localized in different 
wells the A bosons in the left and the B boson in the right well. Similar to the previous 
scenario, increasing qab leads again to an increase of the tunneling period, an effect which is 
intuitive here since the components that are initially prepared in different wells, are forced 
to stay apart from each other by the repulsive interspecies force. Moreover, an important 
point to note is that the contributing states are always those of the lowest band and the 
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Figure 6. (color online) Population of the right well (a) ua of species A and (b) of species B 
at gAB = 0.2 for different values of qa for a species separated initial state (c) Energy spectrum for 
9AB = 0.2 

number states mainly involved are \AA,B) and \B,AA) since the former one is the initial 
state. Therefore as long as qa remains comparatively small, the dynamics consists of a slow 
oscillation between the states \AA,B) and \B,AA) and is correlated in the sense that the 
A and B bosons always occupy different wells in the course of the dynamics. 

The above is shown in Fig. |6] (a),(b) for qab = 0.2, where the population of A and B 
bosons in the right well is plotted. An important difference compared to the completely 
imbalanced preparation, is that the increase of qa leads here to a decrease of the tunneling 
period T initially, reaching a minimum at qa ~ 0.2. Subsequent increase of qa leads to an in- 
crease of the period again. Resorting to the energy spectrum for an explanation (Fig. Et^c)), 
one should focus on the lowest doublets which have dominant contributions from \AA, B) and 
l-B, AA). We see a splitting of the lowest doublet as they approach the avoided crossing lead- 
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ing to an increase of the tunneling rates. For larger gA, it is the energetically excited doublets 
which represent the main contribution. The two levels of the excited doublet come closer 
in energy as increases further leading to a smaller tunneling frequency. In terms of tun- 
neling processes, the dominant sequence here is \AA, B) — ?■ \A,AB) — ?■ \AB,A) — )■ \B,AA) 
and a somewhat suppressed sequence is \AA,B) — \A,AB) — \0,AAB) \B,AA). Us- 
ing Eq. ([3]), the tunneling rates scales as / ~ /2gAB{gAB — 9 a) for the first sequence 
and / ~ J^/{gAB — QaY f*^^ the second sequence. We see that, as gA increases from zero, 
the tunneling frequency increases, reaching a maximum for gAs = gA (it does not actually 
diverge as the formula suggest since in this case other higher order terms with respect to 
J becomes relevant) while beyond this point the period increases again. The crucial dif- 
ference to the previous case of a completely imbalanced initial condition is the sign in the 
denominator which was positive previously and is negative here. This implies that while in 
previous case, there was a monotonic increase in the tunneling period with increasing gA, 
here we encounter an initial decrease as gA approaches gAS- For very high values gA = 25.0 
additional states contribute to the dynamics leading to the high frequency 'noise' observed. 

The avoided crossing present for the lowest lying states of the energy spectrum occurs 
also for higher values of gAB near gA ~ Qab- This results in a similar dynamical behavior 
with respect to the dependence on gA- Different is the case of high interactions gAB = 25.0, 
which is illustrated in Fig. [71 Here the tunneling period decreases substantially as gA takes 
larger values. The very smooth behavior for gA = ^ (Fig. E(a)), where in principle only 
the lowest doublet \AA,B) ± \B,AA) contributes, changes to rapid small oscillations and 
erratic patterns as the intra-species interaction increases indicating that other higher lying 
states are involved in the dynamics. The strong intra-species repulsion here serves again as 
the principal destructor of the correlated shuffling between the initial state and its mirror 
number-state. We can attribute the increase of the tunneling rate, to the increase of the 
density overlaps due to intra- and inter-species strong repulsion, in the line of the arguments 
provided in the discussion of Fig. |H 

As a last remark on the dynamics of the species separated initial state we would like to 
comment on the degree of correlation of the tunneling. Since the tunneling consists here in 
principle of a shuffling between \AA, B) and \B, AA), the two species spent most of the time 
in different wells. Therefore the probability to find B and A species in the same well remains 
always close to zero, while the A particles tunnel as a pair. Similar to the previous section 
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Figure 7. (color online) Population of the right well n/j, ua and ns for qab = 25.0, (a) gj^ = 0.0, 
(b) gA = 5.0, (c) gA = 20.0 for the species-separated initial state. 

this behavior ceases to exist in general for strong qa where single particle tunneling for the A 
species via excited states is induced. Note that for the so-far discussed cases of initial states, 
the destruction of the correlated tunneling behavior, (three bosons staying together; the two 
species remaining separated), results from a strong increase of the intra-species interaction 
which drives the system beyond the simple number state dynamics {\AAB,0) -^^ \0,AAB) 
or \AA,B) \B,AA)). We show next that such strong deviations from the initial state 
configuration can also be achieved for the situation of a partially population imbalanced 
initial state but for a different reason. 



V. PARTIAL POPULATION IMBALANCED INITIAL STATE. 

A novel tunneling mechanism is encountered if the initial state is prepared such that 
the two wells share an equal mean value of the population of A atoms while the B atom 
is on the left well. This initial state we call partially population imbalanced state. The 
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behavior observed above namely the increase of the tunnehng period with increasing qab 
and its decrease with increasing qa can still be observed here. However, a major difference 
compared to the preceding cases arises in terms of the evolution of the different states which 
reflects itself in the corresponding time-evolution of the populations. 

In Fig. [8] we show the populations n^, Ub and n/j for qab = 0.2. Naively, one would 
expect that the the B boson will undergo Rabi-oscillations on the background of the A 
bosons which should remain with equal population in each well. However this does not 
happen for qab > Qa- The envelope behavior of the A particle population i.e ua in Fig. 
IHl^a) for = first increases then decreases, indicating that the single A atom in the right 
well tunnels partially to the left well thus decreasing the population of the A particles in 
the right-well. The B boson on the other hand tunnels completely to the right well. This 
process is retained thereafter and is overall periodic. The envelope behavior is modulated by 
high frequency oscillations of significant amplitude involving a rapid tunneling between the 
two wells. As Qa is increased from to 0.2 the pattern becomes more irregular consisting 
mainly of a constant envelope and shows rapid oscillations. The amplitude of the oscillation 
of ua remains large. When the intra-species interaction strength qa = 0.3 becomes larger 
than the inter-species coupling qab = 0.2 (Fig. El^c)), the tunneling of A bosons is strongly 
suppressed. For even higher interactions qa = 5.0 (Fig. EJ^d)), the A bosons are completely 
localized while the B boson undergoes Rabi oscillations between the two wells as one would 
expect intuitively since the highly repulsive species A are initially in different wells. 

The evolution of the dynamics shows further characteristics for stronger interspecies in- 
teractions. Fig. E] presents the results for qab = 5.0. For qa = 0.0 (Fig. EJ^a)), there are 
two distinct oscillations for both ua and ns'- a fast fluctuation with significant amplitude 
for ub coupled to a large amplitude motion of ua- Intuitively one can understand this 
behavior (seen also in the previous case) for large qab as follows: the tunneling of the B 
boson to the right well pushes the A bosons to the left well due to the strong repulsion and 
vise versa leading to a counterflow type of dynamics. The fast oscillation of considerable 
amplitude for ua involves tunneling of a 'complete' A boson and partial tunneling of a B 
boson between the wells. The origin of these oscillations can be understood via the number 
state decomposition of the initial state as will be explained below. Opposite to this, for 
Qa = 4.0 (Fig. ID^b)), the tunneling of A bosons is considerably suppressed and the B boson 
undergoes a rapid oscillation between the wells. For even higher qa as before we get an 
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Figure 8. (color online) Population in the right well ur, ua and ub for qab = 0.2 for (a) gA = 0.0, 
(b) QA = 0.2, (c) OA = 0.3 (d) QA = 5.0. 

almost complete suppression of the A boson tunneling while the B boson executes the same 
very fast oscillations. 

For very strong inter-species interaction qab = 25.0, a similar pattern is seen for low qa 
(not shown) albeit with a much longer period. For quite strong qa = 5.0 there is a tendency 
for suppression of the tunneling of the A boson (Fig. M^c)) which still oscillates but with a 
small amplitude. Unlike qab = 5.0, increasing the interaction to qa = 20.0 (Fig. M^d)) does 
not reduce the tunneling of the A bosons but increases it approaching a 'fermionization' 
type behavior of the dynamics. 

To identify the underlying dynamical mechanisms leading to the above observations, we 
first note that the initial state in this case is not necessarily a pure number-state \AB, A) but 
is a linear combination of the number states: \AAB,0), \AB,A) and \B,AA) maintaining 
the required population balance of the initial state (equal population of A bosons in each 
well and the B boson in the left). 

For this initial setup the tunneling dynamics consists of transferring the atoms between 
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Figure 9. (color online) Population in the right well n^, ua and ns for (a) g^s = 5.0 and = 0.0, 
(b) gAB = 5.0 and gA = 4.0, (c) gAB = 25.0 and gA = 5.0 (d) gAB = 25.0 and gA = 20.0 . 



the initial state and a target state composed of the number-states \0,AAB), \A,AB) and 
\AA,B). For qab ^ Oa, the number-state \AA,B) represents the dominant contribution to 
the target state and thus the dynamics consists of transferring the atoms between the initial 
state and the configuration \AA,B). As a result we have a tunneling of the B boson to 
the right well and of a single A boson to the left well which we can observe in the envelope 
behavior of ua and of Fig|8](a) and more prominently in Fig. M^a). The faster oscillations 
are the result of the contributions from the states |0, AAB) and \A, AB). For qab ~ gA, we 
have contributions of approximately the same magnitude from almost all the number states 
leading to Josephson like oscillations. However, for qab ^ 9a, the dominant contribution 
of the target state is \A, AB) and thus the system shows a transfer between the initial state 
and the state \A,AB). Therefore the A bosons are effectively localized while the B bosons 
undergo Rabi oscillations between the wells. 
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VI. CONCLUSION AND OUTLOOK 



We have investigated the tunnehng dynamics of a strongly correlated few body bosonic 
binary mixture in a one-dimensional double-well covering the complete range of intra- and 
inter-species interaction. Our focus is the interplay of the inter- and intra-species correlations 
and their impact on the dynamics. We observe that the tunneling period increases drastically 
as the inter-species interaction g^s increases which is due to quasi-degenerate symmetric 
states contributing primarily to the dynamics. This effect is quite general and observed here 
for different initial configurations. 

The intra-species coupling qa possesses a different impact on the behavior of the dynamics, 
depending on the strength qab as well as on the initial state. The general trend is that for 
large qa the overlap of localized wave functions of contributing states becomes larger and 
thus the effective tunneling coupling is increased leading to higher tunneling frequencies. 
For low interactions though different behavior is encountered for different setups. For a 
completely imbalanced initial state, for instance, we observe that for small values of Qab, 
the tunneling period increases as we increase qa in the weak interaction regime. However 
for larger values of qa, the tunneling period reduces with increasing qa- This behavior is 
not seen for the species-separated initial condition. In the latter case, we observe a minimal 
period at qa = Qab which is a manifestation of an avoided crossing in the spectrum. 

Concerning the different initial states of the ensemble the complete population imbalanced 
state exhibits generically a completely correlated tunneling process for the A and B species 
which breaks only for large values of qa leading to an attempted single particle tunneling 
and independent fermion-like behavior. For the species separated scenario the two species 
tend to stay in opposite wells when the inter-species repulsion is large, a behavior which 
alters only if qa becomes also large. For the partially population imbalanced case where the 
mean population of A atoms in each well n^i = 1, although one would intuitively expect 
that the A particles remain in different wells due to their initial preparation, this happens 
only if the interaction between them is considerably large. In the other cases the A particles 
undergo oscillations and the initially mixed state where an A and a B boson coexist in the 
same well can turn into a separated state for which the A and B species reside in different 
wells. 

Understanding the fundamental effects and mechanisms of the tunneling dynamics in 
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strongly correlated bosonic mixtures on a few body level can be seen as a starting point 
to realize systems such as bosonic transistors or to create schemes for selective transport of 
individual bosonic component in reservoir-sink systems as well as for studies of entanglement 
and statistical properties of mixed ensembles. Further considerations could include higher 
number of particles or species and effects of parameters which differentiate the two species. 



Appendix A: Computational Method MCTDH 

Our goal is to study the bosonic quantum dynamics for weak to strong interactions in 
a numerically exact fashion. This is computationally challenging and can be achieved only 
for few atom systeni. Our approach is the Multi- Configuration Time Dependent Hartree 
(MCTDH) method 
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48j being a wave packet dynamical tool known for its outstanding 
efficiency in high dimensional applications. 

The principle idea is to solve the time- dependent Schrodinger equation 



as an initial value problem by expanding the solution in terms of Hartree products $j = 

m = Y.Mmj{t)- (Ai) 

The unknown single particle functions ipj{j = 1, ...,n, where n refers to the total number 
of single particle functions used in the calculation) are in turn represented in a fixed primitive 
basis implemented on a grid. The correct bosonic permutation symmetry is obtained by 
symmetrization of the expansion coefficient Aj. Note that in the above expansion, not only 
are the coefficients Aj time dependent but also the single particle functions ifj. Using the 
Dirac-Frenkel variational principle, one can derive the equations of motion for both Aj and 
$j. Integrating these differential equations of motion gives us the time evolution of the 
system via flAll) . This has the advantage that the basis $j(t) is variationally optimal at 
each time t. Thus it can be kept relatively small, rendering the procedure more efficient. 

Although MCTDH is designed primarily for time dependent problems, it is also possible 
to compute stationary states. For this purpose tKe .etKod,s used Q. The Uey 
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idea is to propagate a wave function \l/o by the non-unitary operator e~^'^ . As r — )> oo, 
this exponentially damps out any contribution but that stemming from the true ground 
state like e~^^'^~-^°^'^ . In practice one relies upon a more sophisticated scheme called the 



improved relaxation l5l| which is much more robust especially for excited states. Here 
(\i/|^r|\E') is minimized with respect to both the coefficients Aj and the orbitals (pj. The 
effective eigenvalue problems thus obtained are then solved iteratively by first solving Aj 
with fixed orbital ipj and then optimizing ipj by propagating them in imaginary time over a 
short period. This cycle is then repeated. 

We note here that the computation of very long tunneling times using the MCTDH 
propagation scheme is numerically impractical. For these cases we computed the dynam- 
ics through the expansion of few-body eigenstates. Moreover, for extremely close quasi- 
degenerate states, convergence is difficult. In these simultaneous relaxation of a 
whole set of these eigenstates keeping them orthogonal is performed by a method known as 
block relaxation. 
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